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Popular summary: 

Aerosols are small particles suspended in the atmosphere and have a variety of natural and man-made 
sources. Knowledge of aerosol optical depth (AOD), which is a measure of the amount of aerosol in the 
atmosphere, and its change over time, is important for multiple reasons. These include climate change, 
air quality (pollution) monitoring, monitoring hazards such as dust storms and volcanic ash, monitoring 
smoke from biomass burning, determining potential energy yields from solar plants, determining 
visibility at sea, estimating fertilisation of oceans and rainforests by transported mineral dust, 
understanding changes in weather brought upon by the interaction of aerosols and clouds, and more. 

This paper evaluates a new dataset which provides AOD over land from Sea-viewing Wide Field-Of- 
View Sensor (SeaWiFS) measurements. The dataset is available free of charge from 
http://disc.gsfc.nasa.gov/dust . The advantage of SeaWiFS is that up to now it provided a longer time 
series (late 1997 to December 2010) of measurements than any other sensor used for determining 
column AOD. The quality of the data is comparable to that from other satellite sensors. 


Abstract: 

This study evaluates a new spectral aerosol optical depth (AOD) dataset derived from Sea-viewing 
Wide Field-of-view Sensor (SeaWiFS) measurements over land. First, the data are validated against 
Aerosol Robotic Network (AERONET) direct-sun AOD measurements, and found to compare well on 
a global basis. If only data with the highest quality flag are used, the correlation is 0.86 5 and 72% of 
matchups fall within an expected absolute uncertainty of 0.05 + 20% (for the wavelength of 550 nm). 
The quality is similar at other wavelengths and stable over the 13-year (1997-2010) mission length. 
Performance tends to be better over vegetated, low-lying terrain with typical AOD of 0.3 or less, such 
as found over much of North America and Eurasia. Performance tends to be poorer for low-AOD 
conditions near backscattering geometries, where SeaWiFS overestimates AOD, or optically-thick 
cases of absorbing aerosol, where SeaWiFS tends to underestimate AOD. Second, the SeaWiFS data 
are compared with midvisible AOD derived from the Moderate Resolution Imaging Spectrometer 
(MODIS) and Multiangle Imaging Spectroradiometer (MISR). All instruments show similar spatial and 
seasonal distributions of AOD, although there are regional and seasonal offsets between them. At 
locations where AERONET data are available, these offsets are largely consistent with the known 
validation characteristics of each dataset. With the results of this study in mind, the SeaWiFS over-land 
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Multiannual average annual cycle of AOD from a variety of different satellite datasets, for different 
regions of the world. This illustrates the differences in AOD in different regions (due to different 
aerosol sources and transport pathways), as well as the general consistency between the satellites’ 
estimates of seasonality. It also, however, highlights remaining uncertainties related to small offsets 
between the datasets; understanding these will help to improve future versions of the data. 
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Abstract. This study evaluates a new spectral aerosol optical depth (AOD) dataset derived from 
Sea- viewing Wide Field-of-view Sensor (SeaWiFS) measurements over land. First, the data are val- 
idated against Aerosol Robotic Network (AERONET) direct-sun AOD measurements, and found to 
compare well on a global basis. If only data with the highest quality flag are used, the correlation 
5 is 0.86 and 72 % of matchups fall within an expected absolute uncertainty of 0.05 + 20 % (for the 
wavelength of 550 nm). The quality is similar at other wavelengths and stable over the 13-year 
(1997-2010) mission length. Performance tends to be better over vegetated, low-lying terrain with 
typical AOD of 0.3 or less, such as found over much of North America and Eurasia. Performance 
tends to be poorer for low- AOD conditions near backscattering geometries, where SeaWiFS overes- 
1 0 timates AOD, or optically-thick cases of absorbing aerosol, where SeaWiFS tends to underestimate 
AOD. Second, the SeaWiFS data are compared with midvisible AOD derived from the Moderate 
Resolution Imaging Spectrometer (MODIS) and Multiangle Imaging Spectroradiometer (MISR). 
All instruments show similar spatial and seasonal distributions of AOD, although there are regional 
and seasonal offsets between them. At locations where AERONET data are available, these offsets 
1 5 are largely consistent with the known validation characteristics of each dataset. With the results of 
this study in mind, the SeaWiFS over-land AOD record should be suitable for quantitative scientific 
use. 
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1 Introduction 


The remote sensing of aerosol optical depth (AOD) from satellites over land is generally more com- 
plicated than over ocean. Reasons for this include a brighter and more heterogeneous surface re- 
flectance, diversity in aerosol microphysical properties, difficulties in cloud identification, and, par- 
ticularly near strong aerosol sources, heterogeneity of the aerosol itself. The limited information 
content of current and previous passive satellite imaging radiometer measurements (Hasekamp and 
Landgraf, 2007) has led to the development of many different AOD retrieval algorithms, utilising 
the strengths and mitigating the weaknesses of each sensor (e.g. Kokhanovsky and de Leeuw, 2009). 

Retrieved AOD from these datasets is typically evaluated using ground-based estimates of AOD, 
from measurements of the extinction of direct solar radiation using sun-photometers (e.g. Ichoku 
et al., 2002, Kahn et al., 2010, Levy et al., 2010, Sayer et al., 2012). The methodology adopted 
in this work draws from these studies. This evaluation is an important component in establishing 
the reliability of the satellite data for scientific use, such as in studies of air quality (Gupta and 
Christopher, 2008, Hoff and Christopher, 2009), radiative forcing (Chylek et al., 2003, Bellouin 
et al., 2005), or trend detection (Zhang and Reid, 2010, Yoon et al., 2011). The difficulty of accurate 
AOD retrieval over land has so far limited some of these applications (e.g. Li et al., 2009). 

At the time of writing, the longest single-sensor record suitable for determination of column AOD 
over land and ocean surfaces is the Sea- viewing Wide Field-of-view Sensor (SeaWiFS), which was in 
operation between September 1997 and December 2010, with a small number of temporary outages. 
The main focus of the SeaWiFS mission was determination of ocean color parameters (McClain 
et al., 2004). Due to the small contribution of the ocean to the reflected solar radiance at top-of- 
atmosphere (TOA), this requires a highly accurate and stable calibration, which was maintained 
over the mission (Eplee Jr. et al., 2011; of order 1-3 % absolute calibration accuracy, and temporal 
stability around 0.3 %). The long time series and quality of calibration make SeaWiFS a good 
candidate sensor for the creation of an AOD dataset. With an equatorial crossing time around noon 
(albeit drifting into the early afternoon in its final few years), SeaWiFS’s sampling time was in 
between other satellite platforms which carry instruments suitable for AOD retrieval (e.g. Terra, 
ERS-2, or Envisat in the mid-morning, and the A-Train in the early afternoon). 

For these reasons, the Deep Blue Utilization of SeaWiFS Through the Data and Information Ser- 
vices Center (DUST-DISC) project has created such an AOD record from SeaWiFS. The current 
version (v003) is freely available from http://disc.gsfc.nasa.gov/dust/ along with supporting docu- 
mentation for users. Both Level 2 (individual orbit) and Level 3 (daily/monthly gridded composite) 
products are available. Further information about the dataset, including details and validation of the 
algorithm used over water, is provided by Sayer et al. (2012). The over-land AOD retrieval draws on 
the heritage of the ‘Deep Blue’ algorithm (Hsu et al., 2004, 2006), which is part of the operational 
MODerate Resolution Imaging Spectroradiometer (MODIS) processing, and will be described in a 
forthcoming publication. The aim of this paper is to evaluate the over-land AOD dataset through 
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validation of Level 2 AOD retrievals with Aerosol Robotic Network (AERONET) data (Holben 
et al., 1998) and comparison with other satellite AOD datasets derived from MODIS and Multiangle 
Imaging Spectroradiometer (MISR) measurements on a Level 3 basis. 

Section 2 describes the satellite and ground-based data used in the study, and Section 3 presents 
some general results concerning the dependence of comparison statistics on factors such as geometry 
and data quality flags. Next, Section 4 discusses individual regions in more detail. Finally, Section 
5 summarises the main results of the study. 

2 Datasets and methodology used 

The over-land SeaWiFS aerosol dataset provides retrieved AOD at the commonly-used reference 
wavelength of 550 nm, at the SeaWiFS bands centred near 490 nm and 670 nm, and, over some 
surface types, the band centred near 412 nm. The primary quantity of interest is the AOD at 550 
nm, 7550 , where the subscripted number indicates the wavelength in nm. Unless otherwise specified, 
references to AOD imply 7550 . Each Level 2 retrieval pixel has an associated quality assurance (QA) 
flag with a value between 1 and 3, with 3 indicating the most confidence in the retrieval and 1 the 
least; over land, only QA=3 retrievals are used to create the Level 3 products. The comparisons 
presented in this paper can be divided into two types, described in the succeeding sections. 

2.1 AERONET validation 

Comparison of spectral AOD with AERONET data on a Level 2 basis helps to establish the reliabil- 
ity of the SeaWiFS data in many different aerosol/surface regimes. The AERONET record provides 
spectral AOD in several bands between the ultraviolet (UV; 340 nm) and shortwave infrared (swIR; 
1,640 nm) with an uncertainty of order 0.01-0.02; the larger uncertainties are found in the UV (Hol- 
ben et al., 1998, Eck et al., 1999). The exact band configuration depends on the specific instrument. 
The validation protocol is as adopted in Sayer et al. (2012), and previously-cited studies. AERONET 
provides a point measurement with a temporal resolution around 15 minutes, while SeaWiFS pro- 
vides a single-time measurement with a spatial footprint of 13.5 km x 13.5 km (at the sub-satellite 
point). Additionally, AERONET samples the direct atmospheric path between the Sun and the site, 
while SeaWiFS measures light which has traveled this path and then been reflected into the view- 
ing direction. To mitigate these sampling differences, the former dataset is averaged in time, and 
the latter in space. A valid matchup is defined when there is at least one AERONET measurement 
within 30 minutes of the satellite overpass, and at least one SeaWiFS retrieval in a 25 km radius of 
the AERONET site. 

Only cloud-screened and quality-assured Level 2.0 AERONET data (Smirnov et al., 2000) are 
used. A total of 183 AERONET sites, which each provide at least 10 matchups when the SeaWiFS 
data are restricted to the highest quality flag (QA=3), are used. Note that only 109 of these sites 
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include a comparison of AOD at 412 nm (due to absence of AERONET and/or SeaWiFS data in the 
relevant bands). 

As the SeaWiFS and AERONET bands do not match, in this study the AERONET AODs are 
interpolated spectrally to the relevant SeaWiFS wavelength using the Angstrom power law, a = 
dlnr/dlnA, where the Angstrom exponent a used for this interpolation is provided in the standard 
AERONET data product, and the interpolation uses the spectrally-closest r and a to the desired 
SeaWiFS wavelength. The spectral variability of AOD is generally smooth, so the error introduced 
by this interpolation is negligible. AERONET data are denoted with a subscripted A, i.e. ta. The 
SeaWiFS dataset provides a evaluated between 490 nm and 670 nm; this is also compared with 
the closest AERONET a (440 nm-675 nm), and unless specified otherwise, further references to a 
indicate this combination. 

2.2 Multi-satellite comparison 

The second type of comparison involves satellite-derived midvisible AOD (550 nm for SeaWiFS and 
MODIS, 558 nm for MISR). In the interests of brevity, and because 550 nm is a common reference, 
only this one quantity is considered. As sampling and aggregation choices can have significant 
effects on such comparisons (Levy et al., 2009, Sayer et al., 2010), to minimise the effects of spatio- 
temporal mismatch between the data, daily Level 3 data are used. The purpose of this is to show 
how the SeaWiFS data compare to the widely-used MODIS and MISR products on a global and 
regional scale. Multiple publications have already addressed the comparison of over-land MODIS 
and MISR data products with AERONET and each other (Liu and Mishchenko, 2008, Kahn et al., 
2009, Mishchenko et al., 2009, Kahn et al., 2010, Levy et al., 2010, Mishchenko et al., 2010, Hyer 
et al., 2011, Kahn et al., 2011, Shi et al., 2011), so this will only be discussed where relevant in 
the context of SeaWiFS. The local overpass times of the satellites vary between late morning and 
early afternoon; Smirnov et al. (2002) observed that variability of AOD between these times is of 
order 10 % of the daily mean AOD or less for most aerosol types, and so is unlikely to be a major 
contributing factor to observed offsets in many regions (although aerosol transport may contribute 
to scatter). 

The comparison uses the standard Level 3 daily products for the period 2006-2010. For MISR, 
this is the version 22 dataset (Martonchik et al., 2009, Kahn et al., 2010). For MODIS, two sets 
of products are considered. Firstly, the MODIS Collection 5.1 standard product, ‘Deep Blue’ over 
arid surfaces (Hsu et al., 2004, 2006) and ‘Dark Target’ over vegetated surfaces (Levy et al., 2007). 
Secondly, the ‘data-assimilation (DA) quality’ MODIS AOD (hereafter DA MODIS), which is an 
attempt to create a filtered and bias-corrected version of the MODIS Collection 5.1 ‘Dark Target’ 
data (Hyer et al., 2011). The MODIS sensors aboard the Terra and Aqua satellites are considered 
separately. 

The comparisons with MODIS are at 1° horizontal resolution, while that with MISR is at 0.5°, 
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1 25 corresponding to the different resolutions of their standard products. Because the mutual overlap 
of all datasets on a given day is low, colocations are performed between SeaWiFS and each other 
dataset on an individual basis. To be considered for the comparison, grid cells in the daily products 
are only considered where both datasets contain at least 10 retrievals for the 1° comparisons, or 3 
for the 0.5° comparisons, these thresholds being roughly equivalent relative to the grid size. 

1 30 2.3 Definition of regions 
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Fig. 1 . Geographical bounds of land regions used in this study for comparison of SeaWiFS AOD with other 
datasets. 

Both sets of analyses consider over-land data globally, and also stratified into one of nine different 
regions, as shown in Figure 1. The choice of regions is somewhat arbitrary, and is not intended to 
follow geopolitical boundaries, but allows a fairly coarse examination of the retrieval under different 
atmospheric and surface conditions. For example, a generally low AOD and dense vegetation are 
1 35 found over much of eastern North America or Eurasia, but an arid surface and periodic dust storms 
are found in northern Africa and the Middle East. Nine regions, corresponding roughly with those 
often used in regional analyses of aerosol properties, were chosen as a compromise between iso- 
lating individual aerosol/surface regimes and retaining sufficient data volume in each to present a 
meaningful analysis. 
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Fig. 2. Scatter density plots between SeaWiFS and AERONET spectral AOD, for (a) 412 nm, (b) 490 nm, (c) 
550 nm, (d) 670 nm, (e) Angstrom exponent (for all points), (f) Angstrom exponent (for ta ,550 > 0.3). Only 
points where QA=3 are included. Correlation coefficients are shown above each panel. 


3 Global comparison with AERONET data 
3.1 Effect of QA 

Considering data from all AERONET sites together, and restricting SeaWiFS data to QA=3 only, 
leads to joint histograms of SeaWiFS and AERONET AOD or Angstrom exponent as shown in 
Figure 2, revealing a high level of agreement in spectral AOD. The Angstrom exponent compares 
poorly, but if only cases of moderate or high AOD (ta, 55 o > 0.3) are considered, the correlation 
improves. This is expected as in low- AOD conditions small uncertainties on r can propagate to 
large uncertainties on a (Wagner and Silva, 2008), and similar behaviour is generally found in other 
satellite datasets (e.g. Levy et al., 2010). Some artefacts are visible in the distributions; namely, that 
the minimum permitted 7550 is set to 0.02 for SeaWiFS for some surface types, and that a is forced 
to 1 in SeaWiFS if 7550 <0.1 (due to a lack of information), and bounds are set on a dependent on 
the aerosol microphysical model used for retrieval in a given area. 

An important concept (Kahn et al., 2010, Levy et al., 2010, Sayer et al., 2012, and others) is 
the notion of expected error (EE); namely, what is the confidence envelope which describes the 
uncertainty on the SeaWiFS AOD? This can be estimated from the comparisons with AERONET, 
assuming that the AERONET data represent an unbiased and comparatively low-error measure of the 
true AOD. The goal is to define some function such that the absolute difference between SeaWiFS 
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and AERONET falls within this EE for one standard deviation (~68 %) of the colocations, and two 
standard deviation (^95 %) within twice this EE. The error ratio (ER) for an individual matchup is 
defined as the difference between SeaWiFS and AERONET AOD, divided by the EE. 

The established over-land EE for MODIS Dark Target and MISR midvisible AODs are 0.05+0. 15 ta , 550, 
and the greater of 0.05 or 0.2ta,55o> respectively. An AOD-dependent component is linked to un- 
certainties associated with aerosol microphysical model assumptions, while an AOD-independent 
component is linked to factors such as uncertainties in surface reflectance modelling. Simulations 
and experience with SeaWiFS data suggest that the EE should be larger, as SeaWiFS is lacking in 
spectral resolution compared to MODIS, angular resolution compared to MISR, and spatial resolu- 
tion compared to both. 

Assessment of the set of matchups with QA=3 suggests that an appropriate set of spectral EE 
are 0.05+0.25ta,412, 0.05+0.20ta,49o, 0.05+0.20ta,55o, and 0.05+0. 15ta, 670 , which are met by 
over 68 % of the points (Table 1). Dependent on wavelength, 90 %-95 % fall within twice the EE, 
indicating there are slightly more outliers than would be expected if these expressions of EE were 
true Gaussian measures of the retrieval uncertainty. Because the AERONET uncertainty is much 
smaller, it can be neglected for the estimation of EE. No EE is estimated for a ; as Figure 2 shows, 
it is qualitatively clustered with the AERONET data for ta, 55 o > 0.3, although is not recommended 
for lower AOD. Additionally, for low aerosol loading, the uncertainty on AERONET AOD can itself 
be significant in calculating a (Wagner and Silva, 2008). 

If the requirement for QA=3 pixels is relaxed, the number of comparisons increases, although the 
correlation and fraction of comparisons within the EE decreases (Table 1). This suggests that the 
SeaWiFS quality flags do provide information about the quality of the retrieval. The true decrease 
in quality with decreasing quality flag may be more extreme than noted here, as a primary reason 
for assignment of a low QA flag is suspection of cloud-contamination. By definition the points 
in the AERONET comparison have passed cloud- screening for both the SeaWiFS and AERONET 
datasets, so the frequency of cloud-contamination in SeaWiFS retrievals is likely higher in the whole 
dataset than the subset with valid AERONET comparisons. For these reasons, only QA=3 data are 
recommended for general use, although QA=2 may also be useful if for a specific application a user 
possesses additional information such as auxiliary cloud masks. Only QA=3 data were used for the 
generation of available over-land SeaWiFS level 3 data products. 

Finally, it is relevant whether or not restriction to QA=3 data introduces a sampling bias into the 
underlying distribution of AOD. Figure 3 shows histograms of the AERONET AOD at 550 nm, for 
the nine regions in Figure 1, for the colocations obtained by each different SeaWiFS QA threshold. 

Each region generally approximates a lognormal distribution (also noted by O’Neill et al., 2000 for 
individual AERONET sites), although peaks and widths differ. The Central/South America region 
has a longer positive tail, which is likely due to seasonal biomass burning aerosol observed at several 
of the sites (notably Abracos Hill, Alta Floresta, and Rio Branco). The histograms of AOD are very 
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(d) Eurasia 
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(f) Southern Africa 



(g) North-east Asia (h) South-east Asia 
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Fig. 3. Histograms of AERONET AOD at 550 nm obtained from SeaWiFS-AERONET matchups. Each subfig- 
ure shows a composite from sites in different regions, as defined in Figure 1. Black lines show histograms for 
matchups with all SeaWiFS QA permitted, red lines those matchups obtained using QA=2 or QA=3 retrievals 
only, and green lines those matchups using QA=3 only. 


similar for all QA thresholds, which suggests that restricting to only QA=3 SeaWiFS data does not 
1 95 cause a sampling bias in the ‘true’ distribution of AOD. The region with the least similarity between 
the QA=3 distribution and others is north-eastern Asia. The majority of the points in this region are 
from urban areas in China, or the Gobi Desert. In these cases it appears that restriction to QA=3 
reduces the relative frequency of the lowest (ta, 55 o < 0.05) and highest (ta, 55 o > 1) AOD. 

3.2 Effect of spatial averaging 

200 A further question is whether the spatial averaging of SeaWiFS data is decreasing the apparent 
error, by averaging out noise in the data. To test this, a more restricted additional comparison was 
performed, using only SeaWiFS retrievals whose footprint included the AERONET sites directly (i.e. 
no spatial averaging) and requiring a tighter temporal coincidence of sampling (averaging within 
±15 minutes rather than 30). The resulting sample size, correlations, and fraction within EE are 
205 shown in Table 1 for each QA threshold: statistics are very similar to the 25 km/30 minute case, 
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Minimum Number of matchups R Fraction within EE 

QA Wavelength, nm Wavelength, nm Wavelength, nm 



412 

490 

550 

670 

412 

490 

550 

670 

412 

490 

550 

670 




25 km, 30 minute spatiotemporal average 





1 

27,555 

67,191 

67,202 

67,075 

0.80 

0.81 

0.81 

0.79 

0.64 

0.62 

0.64 

0.61 

2 

22,942 

56,508 

56,512 

56,512 

0.83 

0.84 

0.83 

0.80 

0.68 

0.67 

0.69 

0.66 

3 

12,823 

34,722 

34,727 

34,728 

0.83 

0.86 

0.86 

0.84 

0.73 

0.71 

0.72 

0.69 



Direct spatial matchup, 15 minute temporal average 





1 

16,908 

44,669 

44,659 

44,503 

0.80 

0.79 

0.79 

0.75 

0.64 

0.62 

0.64 

0.61 

2 

12,629 

32,436 

32,426 

32,427 

0.83 

0.85 

0.84 

0.81 

0.69 

0.68 

0.69 

0.66 

3 

6,224 

17,886 

17,879 

17,881 

0.83 

0.87 

0.86 

0.85 

0.72 

0.71 

0.72 

0.69 


Table 1. Statistics of validation of spectral AOD with AERONET data, for different thresholds on the minimum 
SeaWiFS QA value permitted, and both area-averaged and direct colocations. R is Pearson’s linear correlation 
coefficient. 

indicating that retrieval noise is a minor component, i.e. most errors in a given region are likely 
systematic, linked to characterisation of surface and aerosol microphysical properties at that space 
and time, rather than radiometric. The tighter colocation criteria result in sampling decreasing by 
approximately half. Therefore, as the statistics of the comparison are similar, and the sample size 
21 0 double, the 25 km/30 minute comparisons will be used for the further analysis. 

3.3 Dependence on AOD 

Figure 4 shows the error in the retrieved AOD (r — ta) binned as a function of ta, with lines and 
shading to indicate different percentiles of the error distribution. For low AODs, the SeaWiFS data 
have negligible offset. The distribution of retrieval errors in these cases must be biased positive, as 
215 retrieval of negative AOD is not permitted. As AOD increases above about 0.5, at all wavelengths 
a negative bias of order 15 % is observed. However, for the range of AOD in which the majority of 
data lie (typically 0.05-0.3, e.g. Figure 3), the SeaWiFS data have on average an offset of order 0.02 
or less from AERONET. SeaWiFS reports lower a than AERONET for all AOD; as a is fixed at 1 
for 7550 <0.1 in SeaWiFS, Figure 4 (e) suggests that 1.5-2 would be a more appropriate value in 
220 these conditions, to give a comparatively unbiased estimate. 

One possible explanation for a low bias at high AOD is an underestimation of aerosol absorption. 
However, as TOA reflectance is not a linear function of AOD, heterogeneity of aerosol on a sub- 
retrieval scale could also be a factor, particularly close to strong sources (e.g. Alexandrov et al., 
2004, Shinozuka and Redemann, 2011). Figure 5 shows TOA reflectance stored in the lookup tables 
225 used in the SeaWiFS AOD retrieval (defined as the TOA radiance normalised by the solar irradiance) 
for typical observing conditions at the AERONET site of Banizoumbou in central Africa. In this 
figure, solar and viewing zenith angles were taken to be 32°, and the relative azimuth angle 120°. 
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(b) t 490 


( C ) ^550 



(d) t 670 (e) a 




Fig. 4. Retrieval error (SeaWiFS-AERONET) as a function of AERONET AOD, for AOD at (a) 412 nm, (b) 
490 nm, (c) 550 nm, (d) 670 nm; (e) Angstrom exponent (as a function of AOD at 550 nm). In these, data 
are binned in order of ascending AOD, with a bin size of 500 matchups. In each plot, the red line shows the 
mean error, the solid black line the median, dotted black lines the 16th and 84th percentiles in that bin (i.e. the 
68 % confidence interval), and the shaded grey area the bounds of the 2.5th and 97.5th percentiles (i.e. the 95 % 
confidence interval). Only points where QA=3 are included. 


(a) 412 nm (b) 490 nm (c) 670 nm 





Fig. 5. Response of TOA reflectance to changes in AOD at (a) 412 nm, (b) 490 nm, and (c) 670 nm. Geometry, 
surface, and aerosol microphysical properties were taken as those typically encountered at Banizoumbou by 
SeaWiFS. The red line indicates the TOA reflectance as a function of AOD at that wavelength. The black 
dotted lines show the TOA reflectance for AOD of 0.5, 2, and 3.5. The black solid line shows the average of 
reflectances for AOD=0.5 and AOD=3.5, and the equivalent AOD for this reflectance. 
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Surface reflectances of 0.06, 0.1, and 0.2 at 412 nm, 490 nm, and 670 nm respectively were assumed, 
together with an absorbing aerosol model to represent biomass burning smoke. Suppose an AOD 
retrieval is performed on a pixel split in coverage evenly between a thick plume with and AOD of 
3.5 and an elevated background AOD of 0.5. In this case, the area-weighted average AOD of the 
pixel would be 2. However, because of the nonlinearity, the TOA reflectance of this heterogeneous 
case does not correspond to the TOA reflectance of a homogeneous scene with an AOD of 2, and, 
presented with this reflectance, the algorithm would retrieve an AOD of 1.6, 1.7, and 1.8 at 412 nm, 
490 nm, and 670 nm respectively. It is a philosophical point whether this should be considered an 
error or not, and potentially depends on the user’s desired application of the data, as the satellite 
would retrieve the radiative-average AOD rather than the spatial average. 

Although this is a single example, it serves to illustrate the potential magnitude of such effects. 
At high AOD the increase in reflectance with AOD becomes smaller and the curves begin to flatten; 
the fact that the curvature is of this type suggests that, in a high-AOD case, a negative bias in 
satellite-retrieved AOD would be more likely than a positive bias. This will be an issue for other, 
similar retrieval algorithms, and it could be ameliorated in the future by sensors with higher spatial 
resolution, or assuming some degree of spatial heterogeneity in AOD (rather than homogeneity) in 
the radiative transfer calculations used in such algorithms. 


QA 3, t A 5 50 >0.3 



Fig. 6. Retrieval error ratio (ER) at 550 nm, defined as the ratio of the retrieval error to the expected error, as a 
function of AERONET a. Data are binned in order of ascending a, with a bin size of 500 matchups, and only 
QA=3 points where ta ,550 > 0.3 are considered. The lines and symbols are as indicated in Figure 4. 

As a can be used as a first-order indication of aerosol type (e.g. Eck et al., 1999), Figure 6 
examines the ER at 550 nm as a function of q;a- Only points where ta, 55 o > 0.3 are used, both to 
ensure the utility of a a, and as the notion of a dominant aerosol ‘type’ is more meaningful when 
AOD is elevated significantly above typical background values. There is no clear dependence of 
ER on a: a, suggesting that the underestimate of AOD for conditions of high aerosol loading occurs 
for all aerosol types, although the bounds of the 95 % confidence limit become thinner (i.e. fewer 
extreme outliers) as a a increases (corresponding to a more fine-mode dominated aerosol type). 
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3.4 Geometric dependence 
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Fig. 7. Retrieval error ratio (ER) at 550 nm, defined as the ratio of the retrieval error to the expected error, as 
a function of retrieval geometry. The top row shows points where ta ,550 <0.3, and the bottom row ta ,550 > 
0.3. Plots show dependence on (a,e) the solar zenith angle, (b,f) satellite (viewing) zenith angle, (c,g) relative 
azimuth angle, and (d,h) scattering angle. The bin size is 5° for each angle, and only bins with 5 or more points 
are shown. The lines and symbols are as indicated in Figure 4. Dashed vertical lines show the angular range in 
which the central 68 % of the data lie. Only points where QA=3 are included. 

The ER is shown as a function of retrieval geometry in Figure 7, split into low-AOD conditions 
(ta 5 55 o <0.3) where surface reflectance assumptions are expected to be the dominant source of un- 
255 certainty, and high-AOD conditions (ta, 55 o >0.3), where aerosol microphysical property assump- 
tions are expected to dominate. The distributions of ER show some similarities in both groupings. 
When the atmospheric path length of reflected radiance is small (e.g. high-sun, near-nadir observa- 
tions), the 95 % confidence interval is wider, particularly on the positive-ER side. The contribution 
of the aerosol to the TOA signal will be smaller in these cases, and cloud detection also more dif- 
260 ficult; further, this geometry is most frequently encountered in tropical regions, which more often 
have bright reflecting surfaces (e.g. the Sahara desert). The difference between the mean and median 
ER also becomes larger under these geometric conditions for low AOD, further indicative of stronger 
positive outliers. 

A positive bias in AOD as the viewing geometry approaches backscatter is observable in low-AOD 
265 cases (Figure 7 (d)), but not high-AOD cases (Figure 7 (h)). This suggests that the surface reflectance 
models used in the SeaWiFS retrieval underestimate the ‘hotspot’ effect (e.g. Chen and Cihlar, 1997) 
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observed near the backscattering direction, and is therefore one avenue for improvement in future 
versions of the dataset. However, the bulk of the data lie at smaller scattering angles than those 
where this bias is largest. 

Figure 7 (e) reveals an oscillation of ER with solar zenith angle, for ta, 55 o > 0.3. The reasons 
for this are at present unclear; examining the statistics on a regional basis (not shown) suggests the 
pattern is largely from matchups over the North Africa/Middle East, and South-East Asia regions, 
suggesting it could be linked to errors in dust or mixed dust/smoke phase functions at specific geome- 
tries. Finally, for both low-AOD and high-AOD cases, there is a small tendency for the ER to become 
more negative with increasing viewing zenith angle. However, the most dramatic angle-dependent 
effects appear to be associated with high-Sun conditions and retrievals close to the backscattering 
direction. 

3.5 Temporal dependence 


Sites with long time series, land 



Year 


Fig. 8. Retrieval error ratio (ER) at 550 nm, defined as the ratio of the retrieval error to the expected error, 
as a function of time. Data are binned by year. The lines and symbols are as indicated in Figure 4. Only 
points where QA=3 are included, and only those sites with a long AERONET time series (listed in the text) are 
considered. 

The 0.3 % radiometric stability of SeaWiFS calibration (Eplee Jr. et al., 2011) should translate 
into potential artificial trends in retrieved AOD smaller than 0.01 across the mission. The retrieval 
stability may differ from this if non-radiometric factors (e.g. applicability of assumptions made about 
surface reflectance or aerosol microphysical properties) changes through time. This can be assessed 
by examining the change in ER as a function of time in comparisons against AERONET. This ap- 
proach was previously taken by Levy et al. (2010) when in validating the MODIS ‘Dark Target’ 
dataset. Only the subset of AERONET sites including data spanning at least a ten-year time pe- 
riod (Alta Floresta, Avignon, Banizoumbou, Beijing, Bondville, Bratts Lake, CART Site, Dakar, 
Dalanzadgad, GSFC, HJ Andrews, Ilorin, IMS METU-Erdemli, Ispra, Kanpur, MD Science Center, 
Mongu, Sede Boker, Sevilleta, Skukuza, Solar Village, Wallops) are used. Restricting to sites with a 
long time series helps to ensure that changes in the population of available AERONET sites between 
years do not influence the apparent trend in ER. 
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Figure 8 shows percentiles of the ER distribution as a function of year; there are no obvious trends. 
Performing a linear regression of the mean ER against year gives a gradient of 0.01 12±0. 0074 (in 
units of ER per year); assuming a typical 7550 ~0.15 (the mean and median AERONET values for 
the colocated data are 0.21 and 0.13 respectively) gives an estimate of a total drift of 0.01 14±0.0075 
295 (in units of AOD) over the 13-year length of the dataset. Quoted trend uncertainties are the one- 
sigma uncertainty on the fit. These values are not significant at the 90 % level, suggesting that 
if these long-term sites are representative of the SeaWiFS dataset as a whole, then there is not a 
statistically- significant linear change in the quality of the data with time, which is an important step 
in determining the utility of the dataset for trend analysis. 


300 4 Regional comparisons 


(a) Mean AERONET 550 nm AOD 



(c) Correlation coefficient, 550 nm 



(b) Fraction within expected error, 550 nm 



Fig. 9. Statistics related to AERONET 550 nm AOD comparison at each site considered: (a) mean AERONET 
AOD; (b) fraction of comparisons within the expected error (0.05 + 0.2ta,55o); (c) Pearson’s linear correlation 
coefficient between the matched data; and (d) median (SeaWiFS -AERONET) AOD bias. 

The previous sections have established the utility of the SeaWiFS data on a global basis. However, 
due to differences in aerosol and surface properties (as well as observation geometry) in different 
parts of the world, it is logical to examine the data on a regional basis. This will aid the understanding 
of researchers interested in studies of a regional nature. To this end, Figure 9 maps various statistics 
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Region 

Number 

of 

sites 

Number of matchups 
Wavelength, nm 
412 490 550 670 

R 

Wavelength, nm 
412 490 550 670 

Fraction within EE 
Wavelength, nm 
412 490 550 670 

Eastern North America 

23 

1,557 

5,434 

5,438 

5,438 

0.85 

0.83 

0.80 

0.71 

0.79 

0.83 

0.84 

0.79 

Western North America 

20 

1,011 

4,935 

4,937 

4,938 

0.74 

0.64 

0.60 

0.50 

0.74 

0.68 

0.72 

0.74 

Central/South America 

15 

443 

2,178 

2,179 

2,179 

0.95 

0.96 

0.96 

0.93 

0.75 

0.63 

0.67 

0.73 

Eurasia 

48 

4,413 

8,937 

8,937 

8,937 

0.84 

0.83 

0.82 

0.77 

0.80 

0.82 

0.84 

0.80 

North Africa/Middle East 

27 

3,165 

6,246 

6,246 

6,246 

0.73 

0.82 

0.82 

0.80 

0.60 

0.59 

0.58 

0.49 

Southern Africa 

11 

3 

1,663 

1,663 

1,663 

0.10 

0.81 

0.79 

0.74 

0.33 

0.41 

0.46 

0.52 

North-east Asia 

9 

723 

1,637 

1,637 

1,637 

0.88 

0.87 

0.86 

0.83 

0.79 

0.73 

0.72 

0.66 

South-east Asia 

22 

902 

1,879 

1,879 

1,879 

0.81 

0.81 

0.80 

0.80 

0.71 

0.63 

0.61 

0.53 

Oceania 

8 

606 

1,813 

1,811 

1,811 

0.22 

0.37 

0.35 

0.32 

0.62 

0.70 

0.71 

0.72 


Table 2. Statistics of validation of spectral AOD with AERONET data. Sites are subset according to geographic 
region (Figure 1). R is Pearson’s linear correlation coefficient. 


305 related to the comparison with AERONET data at the different sites. Table 2 presents statistics of 
the SeaWiFS/AERONET comparison on a regional basis. 


Dataset 

Number of grid cells 

R 

Bias 

Corned 

MISR 

1,169,582 

0.76 

-0.027 

0.076 

MODIS (Aqua) 

1,888,154 

0.81 

-0.015 

0.081 

MODIS (Terra) 

1,563,994 

0.80 

-0.019 

0.080 

DA MODIS (Aqua) 

831,335 

0.85 

-0.026 

0.060 

DA MODIS (Terra) 

843,381 

0.82 

-0.029 

0.065 


Table 3. Statistics of comparison between colocated SeaWiFS daily- averaged AOD at 550 nm with that from 
other datasets. Only grid cells over land are considered. R is Pearson’s linear correlation coefficient between 
the grid-box average values. The bias is the median bias between the daily average values, defined such that 
positive values indicate SeaWiFS AODs are larger. The scaled median absolute deviation of the difference 
between the two datasets is denoted <j me d. 

Figure 10 shows the seasonal variability of SeaWiFS 7550 , alongside the difference between Sea- 
WiFS and the other satellite datasets considered. Images from MODIS Terra are omitted for brevity, 
as patterns were found to be very close to those from MODIS Aqua in most cases. The difference 
31 0 images have been resampled to 3° horizontal resolution for ease of visualisation, and offsets between 
the datasets were largely observed to be spatially coherent on these scales. However, all statistics 
presented herein are calculated at full (0.5° or 1°) Level 3 resolution. 

Table 3 shows the correlation, bias, and scaled median absolute deviation (cr me d) about the median 
difference in 7550 , defined cr me d(x) =/3(x — x), where ^ indicates a median quantity and /? is a scal- 
315 ing factor. If the underlying distribution is Gaussian, then cr me d is equivalent to standard deviation 
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Fig. 10. Comparison between seasonal mean colocated SeaWiFS AOD at 550 nm and other satellite datasets. 
From left-right, columns indicate the seasons December- January-February (DJF), March-April-May (MAM), 
June-July-August (JJA), and September-October-November (SON). The top row shows (a-d) SeaWiFS seasonal 
mean AOD. Successive rows show the difference in seasonal mean AOD from colocated data between SeaWiFS 
and (e-h) MISR; (i-1) MODIS (Aqua); and (m-p) DA MODIS (Aqua). 
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for /? = 1.4826, which is assumed here. This quantity is discussed further by Sayer et al. (2012). 
Note that it is not appropriate to consider a fraction within expected error for the inter- satellite com- 
parisons. Reasons for this include the fact that none of the satellite datasets can be considered to 
provide a ‘ground truth’; that some scatter is to be expected because of changes in aerosol loading 
320 between the overpasses of the different instruments; that there are differences in spatial sampling 
between the datasets (even though the methodology aims to minimise these); and that the expected 
error for individual retrievals from each dataset cannot be propagated into a Level 3 expected error, 
as it is uncertain which components of the retrieval error are a result of random error and which 
are systematic biases for any given point. These statistics are provided on a global scale (regional 
325 correlations are provided later), and reveal a high level of agreement, although a relative negative 
offset between 0.015 and 0.029 on a global basis. 
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Fig. 11. Annual cycles of 550 nm AOD from the satellite datasets, for each region. Month names are abbreviated 
using their first letters, in order January-December. 

The annual cycle of 550 nm AOD is shown in Figure 1 1 for each region. In this figure, the MODIS, 
DA MODIS, and MISR data are drawn from the set of colocated daily measurements described in 
Section 2.2. For simplicity, the SeaWiFS data are a simple multiannual (2006-2010) monthly mean 
330 in each region (as otherwise the plots would require one SeaWiFS line to correspond with each of 
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the colocated datasets). Note that this means the spatial sampling may differ between the datasets. In 
particular, DA-MODIS lacks coverage over bright surfaces. This figure also illustrates the similarity 
between Terra and Aqua data, for the MODIS and DA MODIS datasets. 

Finally, Figures 12-20 show joint histograms of SeaWiFS data against the other datasets used 
335 (AERONET, MISR, MODIS, DA-MODIS) for each of the nine regions described in Figure 1. A 
caveat here is, again, that the spatial extent of the area of each region which is sampled by AERONET 
is smaller than that covered by the satellites. The SeaWiFS/AERONET comparison are drawn from 
the set of matchups on an instantaneous basis (25 km/30 minute averaging), while the satellite inter- 
comparisons are from the colocated daily Level 3 data. 

340 4.1 Eastern North America 
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(c) Eastern North America, R=0.73 





(b) Eastern North America, R=0.67 



(d) Eastern North America, R=0.79 



Fig. 12. Scatter density plots between SeaWiFS AOD at 550 nm and other datasets, for the Eastern North 
America region (Figure 1). (a) shows the comparison between SeaWiFS and AERONET for stations in this 
region, (b-d) show the comparison between colocated daily SeaWiFS and MISR, MODIS (Aqua), and DA 
MODIS (Aqua) data respectively. Correlation coefficients are shown above each panel. 

This region is characterised by low AOD (mean ta, 55 o ~ 0.1 at all AERONET sites; Figure 9) 
and largely vegetated surfaces. Despite the low AOD, correlations between SeaWiFS and the other 
datasets are fairly strong, from 0.67 for SeaWiFS/MISR to 0.8 for SeaWiFS/AERONET. The fraction 
of SeaWiFS/AERONET matchups within EE is well in excess of 68 % in this region. All datasets 
345 show a similar seasonality, with spring and summer peaks in AOD and a drop in autumn and winter 
(Figure 11). The offset between SeaWiFS and other datasets is small (within 0.02), spatially fairly 
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coherent (Figure 10) and shows seasonal variation (positive offset in JJA, negative in DJF and SON, 
variable in MAM) such that SeaWiFS exhibits a stronger seasonality in AOD. The consistency be- 
tween sensors of the seasonal variability in the offset may reflect deficiencies in the seasonality of 
350 the surface reflectance assumed in SeaWiFS, as the low AOD means that the dominant uncertainty 
source is likely to be surface reflectance assumptions. For 7550 > 0.4 SeaWiFS underestimates AOD 
relative to AERONET, MODIS, and DA-MODIS, although these cases account for a small minority 
of the data. 

4.2 Western North America 



(c) Western North America, R=0.71 





(b) Western North America, R=0.53 



(d) Western North America, R=0.70 
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Fig. 13. As Figure 12, except for the Western North America region (Figure 1). 


355 The comparisons within this region are similar to Eastern North America, although part of the re- 
gion is covered by bare (arid or mountainous) terrain, over which AOD retrieval is more challenging. 
For the region as a whole, 72 % of SeaWiFS/AERONET matchups are within the EE. The arid parts 
of this region contain AERONET sites where SeaWiFS compares more poorly (notably Frenchman 
Flat and Rogers Dry Lake, both dry lake beds); SeaWiFS also tends to overestimate AOD relative 
360 to AERONET, MISR, and DA-MODIS for these arid areas. The negative offset of SeaWiFS against 
MODIS suggests MODIS may also be overestimating AOD for part of this region, which is consis- 
tent with the analysis of Levy et al. (2010). Seasonality is similar to Eastern North America, although 
here the spring peak in AOD is more pronounced, and is again generally consistent between sensors 
(Figure 11). 
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Fig. 14. As Figure 12, except for the Central/South America region (Figure 1). 


The AOD in the Central/South America region shows regional and seasonal diversity. High AOD 
is observed from seasonal biomass burning in the Amazon (AERONET sites of Rio Branco, Abra- 
cos Hill, and Rio de Janeiro), which are reproduced well in all datasets, leading to a high level of 
correspondence. SeaWiFS underestimates AOD for the highest- AOD cases relative to AERONET, 
MODIS, and DA-MODIS, and overestimates relative to MISR, which is consistent with known val- 
idation results of these datasets (e.g. Kahn et al., 2010, Levy et al., 2010). Direct inter-satellite 
comparisons in the heart of the Amazon are sparse due to frequent cloud cover, with significant vari- 
ability between mid-morning and early afternoon (Meskhidze et al., 2009), which limit sampling. 

Away from biomass burning, lower AOD (ta, 55 o ~ 0.2) is found in urban AERONET sites in 
Mexico and Brazil. Although well-correlated, SeaWiFS has a low offset relative to AERONET for 
these sites. This offset is also observed when comparing with the other satellite datasets, and suggests 
an overestimation of surface reflectance in urban areas in this region. Finally, the southern part of 
this region has lower AOD still (ta, 55 o < 0.1); the correlation between SeaWiFS and AERONET 
is low for these sites, as the range of AOD encountered is comparable to the uncertainty on the 
retrieval, although the fraction within EE is large. The area around Buenos Aires is an exception to 
this, as SeaWiFS has a positive offset relative to the other datasets. 
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(d) Eurasia, R=0.82 



Fig. 15. As Figure 12, except for the Eurasia region (Figure 1). 


4.4 Eurasia 

Similar to Eastern North America, AERONET sites within Eurasia are largely found in suburban or 
agricultural regions with a low-AOD background, vegetated surfaces, and AOD peaks in spring and 
385 summer. Again, 84 % of matches with AERONET are within the EE. The offset between SeaWiFS 
and other datasets again varies seasonally within the region, generally positive for the more vegetated 
parts and negative for the less vegetated parts, but is typically of order 0.02 or less. It is most 
pronounced in springtime, where SeaWiFS AOD is higher than the other satellites by up to 0.05, 
but otherwise the seasonality of AOD tracks well. This region contains almost 25 % of all global 
390 matches between SeaWiFS and AERONET, so weighs heavily in the global statistics. 

4.5 North Africa/Middle East 

The aerosol burden in this region is largely composed of wind-blown mineral dust (with particular 
hotspots in the Bodele Depression and western Sahara Desert), with additional contributions from 
biomass burning and anthropogenic emissions (Prospero et al., 2002, Knippertz and Todd, 2012). 
395 The majority of the surface is bare and bright, although the Sahel and savannas towards the southern 
end are more vegetated, as are some coastal portions. It should be noted that coverage of DA- 
MODIS is particularly sparse in this region, as the DA-MODIS dataset is drawn only from the 
‘Dark Target’ MODIS product, which lacks coverage over bright surfaces (as noted in Section 2, 
‘Deep Blue’ data provide coverage in arid regions in the standard MODIS product). Retrieval of 
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Fig. 16. As Figure 12, except for the North Africa/Middle East region (Figure 1). 


AOD from satellite data is challenging in this region due to the bright surfaces, mountainous terrain, 
complicated aerosol microphysical properties (e.g. Kalashnikova et al., 2005, Dubovik et al., 2006), 
and remoteness of locations meaning that many parts of the region are not well-covered by ground 
measurements (Figure 9). Because of the sampling issues, interpretation of the seasonal cycle in this 
region is difficult, although SeaWiFS, MISR, and MODIS all show peaks from May to July and a 
minimum from October to February, with a range of around 0.2 through the year. 

Although well-correlated overall, at several of sites the proportion of SeaWiFS/AERONET matches 
within EE at 550 nm is 0.5 or less, in which cases SeaWiFS underestimates AOD. For the region 
as a whole, it is 0.58. Spatial patterns of the relative offset between SeaWiFS and other satellites 
vary seasonally, up to ±0.3 in some locations, and are not always consistent between the sensors. 
In general, for the highest AODs, SeaWiFS is higher than MISR, similar to MODIS, and lower 
than DA-MODIS, albeit with significant scatter in all cases. Examination of the surface reflectances 
retrieved by MISR and assumed by SeaWiFS in this region (not shown) reveals spatial patterns in 
close correspondence, which suggests that this is not the primary cause for regionally- varying offsets 
between the two. SeaWiFS AOD is lower than MISR over mountainous terrain; retrieval algorithm 
uncertainties related to elevation are likely to affect SeaWiFS more strongly than MISR (as the for- 
mer relies more heavily on violet and blue wavelengths), so this may indicate a low-bias of SeaWiFS 
over mountains. SeaWiFS being higher than MISR for high AOD is consistent with known MISR 
validation results (Kahn et al., 2010). Comparing SeaWiFS and MODIS over the desert is less in- 
structive because they share the ‘Deep Blue’ algorithm approach (although the recent application to 
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SeaWiFS introduces some updates over the algorithm applied in the Collection 5 MODIS product, 
particularly in terms of surface reflectance). 

Due to the difficulty in accurate AOD retrieval here and sparseness of ground measurements, it 
is difficult to assess with confidence which dataset(s) are closest to the truth, and it is likely that 
without more advanced sensors and/or algorithms this will continue for the near future. Users are 
therefore advised to be aware of the differences between satellite AOD datasets in this region. 

4.6 Southern Africa 
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Fig. 17. As Figure 12, except for the Southern Africa region (Figure 1). 


Southern Africa is more influenced by biomass burning and less by mineral dust than the North 
Africa/Middle East region (Swap et al., 1996, Prospero et al., 2002, and references therein), although 
similar considerations about semiarid and elevated terrain, and complexity of aerosol microphysical 
properties, remain. Only 46 % of SeaWiFS/AERONET matches at 550 nm are within the EE. The 
most poorly-performing sites in this region are biomass burning sites in Zambia, where SeaWiFS 
underestimates AOD by 0. 1-0.2; it is possible that some of this is due to heterogeneity of aerosol 
near sources (e.g. Section 3). However, Figures 10 and 17 reveal that aside from at low-lying coastal 
regions, SeaWiFS AOD is generally negatively offset by between 0.05 and 0.25 compared to other 
satellites and AERONET even in low- AOD cases, which suggests an overestimation of surface re- 
flectance. This could also indicate that the simple way in which terrain elevation is treated in the 
current version of the algorithm is inadequate in this region, where the majority of the terrain is at 
an altitude of 1 km or more above sea level. The same comments as made for seasonality of AOD in 
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Central/South America are applicable here. 


4.7 North-Eastern Asia 
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Fig. 18. As Figure 12, except for the North-Eastern Asia region (Figure 1). 


AERONET sites in this region are split between fairly high- AOD regions in industrialised China, 
the Gobi Desert, and vegetation in Siberia; coverage is sparse for these latter two cases. On the 
whole, 72 % of matches are within the EE. Performance at Irkutsk in Siberia is poor (50 % in EE), 
but this is in a valley in mountainous terrain, and only 10 matchups were obtained. At Yulin, only 
28 % are within the EE; this is an elevated site, within a city surrounded by arid terrain. Again, 
at these sites, SeaWiFS is more likely to underestimate than overestimate AOD. The majority of 
points are from sites where the AOD is dominated by urban pollution and dust, at which the level of 
agreement is good, contributing to the correlation with AERONET of 0.86 for the region as a whole. 

Due to seasonal cloud, snow cover, mountainous terrain, and bright desert surfaces, the majority 
of satellite data points in this region are in northern China and North/South Korea, with much of 
Mongolia, Japan, and Russia more sparsely covered. There is little net bias between SeaWiFS and 
MISR or DA-MODIS, although spatially-coherent biases in parts of the region contribute to scatter 
in Figure 18, particularly in the comparison with MISR. SeaWiFS has a more negative bias relative 
to MODIS, with the strongest outliers in the Taklimakan Desert. SeaWiFS exhibits a similar, but 
less negative, offset relative to MISR in this location. Unfortunately, there are no AERONET data 
available in the Taklimakan. SeaWiFS, MISR, and DA-MODIS all show a peak in AOD from March 
to May around 0.25, followed by AOD around 0.15 later in the year. This spring peak is strongly 
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enhanced in the MODIS data. MODIS and DA MODIS also have a secondary peak in December 
and January, absent from the other satellite datasets. 

4.8 South-Eastern Asia 
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Fig. 19. As Figure 12, except for the South-Eastern Asia region (Figure 1). 


Much of this region consists of agricultural land, interspersed with densely-populated cities, in 
which most of the AERONET sites are found. The aerosol loading is linked with large-scale mete- 
orology, and seasonality differs between the Indian subcontinent and the eastern end of the region. 
In the Indian pre-monsoon (April-June), westerly winds blow dust across the Indo-Gangetic plain, 
where it is trapped by the Himalayas; the winter monsoon is associated with fog and thick haze 
from urban pollution and biomass burning aerosols (Prospero et al., 2002, Gautam et al., 2007, 
2010, 2011). Over the Indo-China peninsula, intense biomass burning during the pre-monsoon 
(February- April), in combination with urban emissions, leads to an optically-thick absorbing haze 
layer (Carmichael et al., 2003, See et al., 2006, Gautam et al., 2012); monsoon rains during the 
summer washout much of the aerosol, and cloud cover limits SeaWiFS’s coverage. Similarly, near- 
persistent cloud cover means SeaWiFS coverage over Indonesia is sparse. 

SeaWiFS is well-correlated with the other datasets in this region (Figures 9, 19), although only 
61 % of matches with AERONET are within the EE, as SeaWiFS has a low relative bias, for all 
AOD. Compared to the other satellites, SeaWiFS tends to retrieve lower AOD at those times and 
locations when very absorbing aerosols are present, which, in combination with the AERONET 
results, suggests the aerosol microphysical models used may be insufficiently absorbing. SeaWiFS 
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retrieves higher AOD than MISR over much of the Indo-Gangetic Plain in all seasons; the low bias 
of MISR with respect to AERONET and MODIS over parts of this region has previously been noted 
(Kahn et al., 2009, 2010). Over the region as a whole, SeaWiFS, MISR, and DA MODIS show 
480 AOD higher by 0. 1-0.2 between March and August, relative to the September-February period. As 
in North-Eastern Asia, this peak is much stronger in the MODIS data. 

4.9 Oceania 
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Fig. 20. As Figure 12, except for the Oceania region (Figure 1). 


Although the correlation between SeaWiFS and other datasets is low in Oceania (Figures 9, 20), 
this is because the AOD is reported as persistently low by all datasets, such that the typical range 
485 of AOD encountered is similar to the typical spread of systematic biases over the region. AOD 
varies between around 0.04-0.05 (March- August) to 0.06-0.12 (September-February); offsets be- 
tween datasets are larger during this second period. Overall, 71 % of SeaWiFS/AERONET matches 
are within the EE. The one poorly-performing site is Birdsville, at which SeaWiFS overestimates 
AOD during all seasons, suggesting an underestimate of surface reflectance. Comparisons with 
490 MODIS and MISR reveal a negative offset of SeaWiFS in DJF and SON over the desert regions 
central and western Australia (of order -0.1); however, in comparison with the DA-MODIS dataset, 
the offset is near-zero in many locations, and data are absent in others. Unfortunately, there are no 
AERONET sites in these remote locations. 
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5 Conclusions 


As the remote sensing of aerosol optical depth from space continues to be a challenging problem, it 
is useful that new datasets be validated against ground-truth data, and compared with other similar 
satellite datasets, to gain an understanding of their strengths and shortcomings. This study has 
performed such an evaluation for a new 13-year AOD dataset from SeaWiFS over land (the over- 
water component of the dataset is discussed in Sayer et al., 2012). The main results of the study are 
as follows: 

- On a global basis, 72 % of comparisons between SeaWiFS and AERONET AOD at 550 nm 
had an absolute difference within the expected error of 0.05 + 0 . 2 r^ ?5 5 o, when only data with 
the highest SeaWiFS quality flag (QA=3) are considered. The analysis suggests that restricting 
to QA=3 does not significantly affect the way the underlying distribution of AOD is sampled, 
and that there is not a significant change in the quality of the SeaWiFS data through the mis- 
sion. Performance is similar for the other wavelengths at which AOD is retrieved. 

- Performance tends to be better over vegetated, low-lying terrain with low and moderate AOD, 
such as found over much of North America and Eurasia. Performance tends to be worse 
for low- AOD conditions near backscattering geometries (notably for scattering angles greater 
than 170°), where SeaWiFS overestimates AOD. 

- SeaWiFS tends to underestimate AOD for optically-thick cases of absorbing aerosol, although 
often this underestimate remains within the expected error. 

- On a regional basis, the satellites tend to produce similar seasonality to each other, and to 
AERONET. Differences between the satellite data are generally consistent with the known 
characteristics of the datasets. Unfortunately, many of the regions where large offsets exist 
between the satellite datasets are remote and not well-sampled by AERONET (also noted by 
Shi et al., 2011). 

In conclusion, the SeaWiFS record provides a useful complement to other satellite AOD datasets 
and is suitable for quantitative scientific analysis. It is recommended that users consider only those 
retrievals identified as QA=3 (which is automatically the case with the Level 3 products), in order to 
extract the most reliable data for analysis, and bear in mind the results of this study depending on the 
geographic region of the user’s desired application. Given the high quality of SeaWiFS calibration, it 
is hoped that through further algorithm development some assumptions related to surface reflectance 
and aerosol microphysical properties can be refined, improving the utility of future versions of the 
data. 
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